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ABSTRACT 

We study transient spiral structures in an isothermal, thin, galactic disc. We 
find no instability that can grow everywhere from infinitesimal disturbances, but 
spiral structure does grow in the disc due to an arbitrarily strong, asymmetric, 
central source. An initially finite spiral undergoes transient swing amplification 
as it is gradually wound-up by differential rotation. An independent sequence 
in negative time describes a leading spiral swinging to a trailing spiral. The 
dynamical coupling is established between the swinging potential and the arm 
particles, by ensuring that this potential constrains a locally rotating distribution 
function centred on the arms. This swing amplification propagates in radius at 
the constant rotational speed of the disc, and leaves behind an exponential density 
decline in space and time. 

Subject headings: galaxies-spiral, galaxies-structure, gravitation, spiral arms 



-3 - 



1. Introduction 



In this paper we study logarithmic spiral structure in an isothermal or 'Mestel' disc. 
The axi-symmetric disc is necessarily present as is the axi-symmetric halo, but we ignore 
the non-axially symmetric halo for clarity. It generally has a higher order effect on the disc 
structure. Non-axially symmetric halo structure, consisten t with the d i sc spi ral structure, 



has been studied in a paper that is available on the Arxiv ( Henriksen I ( 2012h ) . 



The spiral structure that we describe in this paper propagates in radius and 
azimuth, and so might be described as a spiral 'wa ve'. However it is not a propa gating 
wave in the sense o f the w ell known linear theory (jLin. C.C. fc Shu. F.H. 



1966): 



sec 



Binney fcTremaine I (120081 1 for a description of later developments). Rather, the non-axially 
symmetric collisionle ss disc particles are co-movin g i n the mean with the spiral potenti al. 
Recent simulation s, (Wada.Baba fc Saitoh"! ( 201lh ). ( iKawata.Grand &: CropperT f 201 lh ) . 



and observations (iFoyle et al 



( 201lh ) encourage this point of view. It is this co-motion of 
the particles w ith the moving arm that leads us to infer an analogy to the classical 'swing 



amplification' ( iToomrel f 1982h ). 



The evolution of the finite spiral potential is in terms of a variable ( = Vt/r, where 
V is the constant disc velocity and t, r have their usual temporal and radial meanings. 
Consequently recurrent structure in ( generates recurrent structure in time and space. 
We find that a given initial spiral evolves through a steep transient density amplification 
near a particular value of (, as it swings and winds up. Subsequently this amplified 
structure dissipates exponentially as ( increases. Consequently the recurrence of similar 
structure persists only for a finite interval or range of radius. Afterwards the process 
must be re-started either by external disturbance or by some internal linear instability 
('self-excitation') if the disc is isolated. 



Such self-excitation behaviour is well known from extensive simulations (e.g. ISellwoocTI 
(1201 lh ) of iso l ated d iscs, even in the t apered Mestel disc that is predicted to be linearly 
stable fjZang I fjl976h ). fjEvansfcRead I fjl998bh ) except for m = 1. However the full Mestel 
disc, which is the self-sim ilar case, is not expected to have well-posed normal modes 
( GoodmanfeEvansl (jl999l )) in the absence of inner and outer boundaries. It is natural then 
to expect the full Mestel disc dynamics to depend only on the self-similar variable £, 
this reflects locally the absence of boundaries. 



since 



In agreement with prior expectations we do not find a truly self-excited instability 
originating from infinitesimal perturbations. There is instead an indefinitely growing spiral 
structure in the disc in terms of the variable ( that originates from the singular centre of the 
system. We infer this because although the disturbance is zero everywhere at t = = £, it 
becomes infinitely strong at the origin for any t > as ( — >■ oo. This singular centre should 
be replaced by a realistic central spiral in a centrally tapered Mestel disc in order to study 
it as the origin of disc spiral instability. This calculation has not been done in this work. 

Despite this central singularity there may nevertheless be an indication of a genuine 
instability in the full Mestel disc, if the spiral wave may be regarded as reflected at large 
radius (where ( = 0). This is implied by the existence of a complementary solution that 
can be interpreted as due to spirals that have passed from trailing to leading after passage 
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through the centre. However in this paper we do not pursue this possibility, and rather 
focus on the evolution of an initially given spiral arm. 

We may understand an initially given spiral arm as the ultimate development either of 
an externally produced disturbance or of some self-excitation. That is, the arm is comprised 
of particles that have been entrained by the wave and subsequently move collectively with 
it. As remarked above such an e mergent wave ha s been found to arise from an unexpected 
instability in recent simulations ( ISellwoocTI ( 2012[ )) of an isolated Mestel disc. It is suggested 



in that work that the instability arises from reflection of the wave at its inner Lindblad 
resonance, where the distribution function (DF) is modified. Our self-similar structure traps 
the particles into a narrow velocity distribution about local co-motion, so that the linear 
Lindblad resonances do not exist. These are replaced by the extended particle trapping 
in the winding non-linear potential, which is the reason for the swing amplification. Thus 
there is no 'cavity' essential to the amplification studied here. 

The DF of the trapped particles is found from the Collisionless Boltzmann Equation 
(CBE), as written in the locally co-rotating frame established by the axi-symmetric disc 
and halo. It is required to share the isothermal self-similarity of the Mestel disc, and it is 
necessarily time dependent due to the differential rotation. The trapping potential must 
also take a specific isothermal self-similar form. 

The details of the preceding assertions are the subject of the following sections. In 
the next section we introduce the formalism for a compatible co-moving, non-axially 
symmetric isothermal DF and potential. This structure is freely falling in an axi-symmetric 
background disc and halo. In a following section the coupling between the spiral potential 
and the comoving DF of the trapped particles is deduced, and the potential solution is 
found. Finally the spiral wave evolution is analyzed graphically and analytically. 



2. Transient, Corotating, Spiral Structure in the Disc 

In this section we construct spiral potential/surface density in a background isothermal 
disc. Such a disc rotates with a constant mean circular speed = V. We know that this 
can not be a steady configuration because non axially symmetric structure winds up in time 
due to the differential rotation. We treat the subsequent time dependence explicitly. We 
find first a distribution function (DF) for the collisionless particles that are comoving with 
the spiral arms and then give a key equation for the compatible potential. 

We treat the spiral structure by remaining close to an isothermal self-similar evolution 
in time, at least before major winding has occurred. This is consistent with the initial 
self-similarity imposed by the Mestel disc. For the collisionless particles the explicit 
collisionless Boltzmann equation (CBE) equation is written in a differentially rotating frame 
with angular velocity Q = V/r. This places our analysis in a reference frame that is in free 
fall in the axi-symme t ric gr avitational background due to the disc and halo. The equation 
becomes dHenriksen I ( j2012h ) 
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d t F +v r d r F + - tv r d r n)d6F 
r 



+ + 2^ + Q 2 r - d r $ dr )d Vr F 



r 



- + 2Qv r + v r rd r Q + -d^ dr )d V4> F = 0, (1) 

where F(r, <p, v r , v^; t) is the two dimensional DF. Other quantities have their usual meaning 
except Q dr , which designates the non-axially symmetric potential of the disc, labelled 
'disc-rotating'. 



The formal self-similar procedure has been discussed elsewhere ( )Le Delliou. Henriksen fc MacMillan 



( 1201 it ) . and references therein) so we will only outline it here. The self-similarity is con- 
strained to be isothermal by the existence of the constant rotational velocity V. We use a 
logarithmic time T as the self-similar Lie parameter and introduce on dimensional grounds 
the scaled quantities R, Y, £, ^ and P according to 

at = e aT , r = Re aT , £ = <p + eT, a = Se" aT 

F = P(R, e, v r , vf, T)e~ aT , v = Y , 

$ dr = ^ In (aR/V) + $$6T + * dr (R, (, 9; T). (2) 

Formally a and e have the dimension of reciprocal time, but in fact all temporal and spatial 
quantities (and consequently velocities) may be thought of as numerical values in terms of 
some fiducial radius r Q and fiducial time t Q when convenient. We have written the potential 
explicitly in terms of these variables, but by using the various definitions its form is seen to 
be equivalent to 

$ rfr = ^ In (or/10 + *dr(R, £, 0; T). (3) 
This is the most general form of an isothermally self-similar, non-axially symmetric 
potential. The coefficient $^ is a constant while the function fy dr contains the angle, 
time and radial dependence. The potential corresponds to the unbalanced non-axially 
symmetric potential in the locally rotating or 'freely-falling' frame. It must satisfy the 
Laplace equation outside the disc and be consistent with the surface density on the disc. 

Through these transformations the independent variables t, r, <p are replaced by the 
scaled variables T, R, £, while the surface density a is replaced by S and the scaled DF 
becomes P. The replacement of v by unsealed Y is peculiar to the isothermality, and is 
used only to emphasize that the velocity is part of the self-similar scheme. The advantage 
of using these variables is that, unless the scaled DF P and pot ential ty dr are strictly 



i ndepe ndent of T, the system is general rather than self-similar (ICarter fc Henriksen 



( 1l99ll )). This allows us to start with self-similar structure and to follow its evolution as the 
time dependence arises. This time dependence arises from the winding of the structure. 

The CBE must be written entirely in terms of these variables and it becomes 
successfully 

aP = d T P+{Y R -aR)d R P + {t+^ + ^§)d i P 
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- ^ + Y w + k d ^ )dY ^ (4) 

with no loss of generality. The explicit dependence on T must vanish for isothermal 
self-similarity. 

The solution for P follows from the characteristic equations 

dP dR 

—— = aP, —— = Y R — aR, 

dT ' dT 

dYn Y* 2VY+ V* 

It - ~R + — + p7-ir- dR * dr > 

dY 1 

= -^(Y,Y R + VY R + d^ dr ), 

where we have written the potential $ dr explicitly. 

We observe that if t, T <— —t, —T and a <— —a, the linear to logarithmic time 
transformation remains formally unchanged together with the transformations to scaled 
variables. This is not merely the time-reversed solution at positive time because of the 
change in sign of a. The origin of the solution at t = is now at T = +oo and it unfolds in 
negative time as T — >■ — oo. Although the characteristic equations are formally the same, 
they describe a different development because the velocities and the quantity e are not 
changed in sign. They would be changed under a strict time reversal. 

The Y^ characteristic equation may be combined with the characteristic expression for 
dR/dT to give 

s( h <« + ^))=-^«. (6) 
which in physical variables is the angular momentum equation 

j t (r(v^ + V)) = -d^ = -d^ dr . (7) 

Moreover the R, £ characteristics may be combined with the Y R , Y^ characteristics to obtain 
an energy equation in the form 

dE' /Q , T . „ Q , T . , Q . T . , , (V\ Y R 



dr 

dT 



(d T V dr - aRd R V dr + ed^ dr ) + (J^J ^d^ dr (8) 

+ ViV + v^^kiRe^y 

Here E' dr = Y 2 /2 + $ dr , where as before $ dr = InaR/V + a^T + ^ dr , is the energy 
in the locally rotating frame. 
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We may also eliminate Yr/ R = d(lnRe aT )/dT between this energy equation and 
equation (jSJ) to obtain 



dT 



(9) 



We note that E dr = E' dr + V(V + Y^), which is the energy equal to $ dr + (Y^ + V r ) 2 /2 + Y%/2 
in the inertial frame but for a constant — V 2 /2. Thus the last equation can be written as 
dEd r /dT equal to all the terms on the right that involve ^^r- 

A swing amplifier exists between the spiral arm and the particles, if the distribution 
function is a function of E' dr . For in such a case the particles are constrained to remain 
near the arm by the potential $d r , and so co-move with the arm on average. We see 
from equation ([9]) that we approach this result by requiring all the potential terms on the 
right hand side to vanish. Once this is applied to equation the inertial energy E& r 
becomes an exact integral, so that the DF is steady in the inertial frame. In order to have 
E' dr as an approximate integral in the locally rotating frame, we require in addition from 
E^ = E' dr + V(V + Y<fr) that <C V. This is normally the case for the majority of particles 
in spiral galaxies. 

This condition on the potential 

d T ^ dr - aRd R V dr + (e + ^-i^§ - a)) d^ dr = 0, (10) 

an n 

is a linear partial differential equation in the potential that imposes a general form . 
The apparent dependence on Yr in this equation is only possible for the potential when 
we remember that the equation holds along a characteristic of the CBE. Hence we may 
eliminate Yr using Yr = aR + dR/dT from the characteristic equations fl5]). The partial 
derivatives are evaluated on the CBE characteristic after having been evaluated holding the 
hidden pair of (T,£,R) constant. 

The equation so written may be solved by using its own characteristics to find that 
^dr = ^dr( r , c) where c = £ — eT + VJ (aR). However from our previous definitions of £, R 
and T one sees that c = <ft + Vt/r = <fi + Q(r)t, which is just the inertial azimuthal angle. 
Thus, as it must be, the potential is steady in the inertial frame where the energy is a 
constant of the motion. 

To obtain a potential that is at least initially compatible with evolving spirals we take 
it to depend on a spiral combination k of c and r plus r as 

= 0- (ii) 

where 

e , e V 
k = cH — In r = 4 H — in R H , 

a a aR 

= <f> + (-)1nr + n(r)t, (12) 
a 
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and note that ( = Q(r)t. 

That this is a solution may be confirmed by calculating (using R) as given in 
equation f lT2|) ) 

(d T ^ dr ) R £ = ard r ^ dr} 

-aRd R ^ dr = -aR{ — — + -^)<9 K * rfr - aRd r V dr , (13) 

aR z aR 

(«+^)*** = 

The sum of these terms is is (V /(a R 2 ) (dR/dT) + aR)d K ^ dr , which is zero along along the 
potential equation characteristic where dR/dT = —aR. Along the characteristic of the 
CBE however dR/dT = Y R — aR, so that to be consistent with the CBE we must regard 
Y R /aR as small. 

In fact a more direct approach is simply to neglect Y R / (aR) compared to one in equation 
( TTOl) . Using the definitions of R and T, this is equivalent to setting (Vt/r)(Y R /V) < 1. As 
the rotating frame radial velocity is small compared to V, this should be well satisfied for 
sufficiently small ( = Vt/r = V/(aR). It transpires that either approach leads ultimately 
to the same form for ^> dr with the latter approximation giving ^ dr (r, k) directly. 

In all this discussion we have supposed for brevity that radii are in terms of a fiducial 
quantity r a = V/a, and we recall that Q = V/r. The second independent variable has the 
form r = Re ST in terms of scaled variables. It is through a dependence on this variable 
that the explicit dependence on T enters, which implies the gradual breaking of the initial 
spiral self-similarity. The equivalence ( = Q(r)t = V/(aR) is helpful when interpreting the 
behaviour found below. 

With this form of the potential and the particle velocities in the locally rotating 
frame being small, equation ([9]) gives E' dr as an approximate integral of the CBE in the 
locally rotating frame. Since E' dr is an approximate integral we expect from the first of 
the CBE characteristics © that P = F(E' dr )e aT . But E = / P dY R dY^ and this should 
be independent of T for initial self-similarity. Thus, recalling the form of the potential in 

scaled variables as $ dr = $^ \n(aR/V) + ct^T + $? dr plus E' dr = Y 2 /2 + $ dn we see that 
an isothermal DF in the locally rotating frame such as 

F(E' dr )=K dr exp[-^), (14) 



removes the T dependence until it develops in the asymmetric potential $ 



dr ■ 



Here K dr is the normalization for the transient, asymmetric, collisionless distribution 
function. The mean velocity of these particles is zero in the locally rotating frame, due to 
the symmetry of the isothermal DF. For the approximation \Y\ < V to hold we require 



(15) 
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To the extent that this rather weak condition holds, the isothermal DF should continue 
to be a solution of the CBE in the potential given by equation f fl8|) below. By treating the 
DF of the particles we avoid detailed consideration of their orbits, but the characteristic 
equations (jSJ) do provide these if necessary 

The potential form (jlip must be made compatible with the Laplace equation above 
the disc, namely 

l(d r (r 2 «9 r ($)) 

+ _l^(sin^$) + -i-^$)=0. (16) 
sin 9 sin 9 

We do this by assuming a 9 dependence (9 is another scaled variable since it is dimensionless) 
in ^dr- After inserting the alternate form for the asymmetric potential (the axi-symmetric 
disc potential already satisfies the Laplace equation above the disc) namely 

$ dr = liir + **.(«, 0,r), (17) 
the Laplace equation becomes the master equation for our purposes in this paper 

+ -d*** + d r (r 2 d r V dr ) 
a 

+ (- )rd r d K ^ dr 

a r 



vt 

a r ' sm 2, 9 
' a(sinaf dr )=0 



dr 



sin 9 



In the Laplace equation the variables t and r are independent. Hence we can replace 
r everywhere by ( = Vt/r since rd r = — Qd^ and d r (r r d r ) = ( 2 d 2 . we will use this form in 
finding the solution for ^dr- We can also note in passing from this equation that one can 
neglect the r dependence in ^ dr at any given t, and so preserve strict initial spiral symmetry 
(k dependence only), only if ( = Vt/r = Q(r)t < e/S. This is generally as expected, but the 
dependence on the tangent of the initial winding angle e/5 is of interest. 

The potential equation ( Fl8|) is linear and can be solved in terms of base functions in 
the form 

= ^lnsm9 + e^T m {9)n m {C). (19) 

We recall that the phase k = + (e/S) lnr + Vt/r and is in the locally rotating frame. 
Naively, by ignoring the ( dependence in the potential, the log spiral is completely wound 
up at a fixed r when Q(r)t = 2tt. This gives ~ 10 7 ' 5 years at r = 10 kpc and V = 200 
km/sec. However even in this naive limit there is an outward moving 'winding wave', given 
by C = Vt/r = est (e/5) lnr, outside of which the pure log spiral remains recognizable 
and similarity is maintained. But as time progresses the ( dependence can not be ignored 
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and it is the consequence of this evolution away from the initial self-similar spiral arm that 
will occupy us in this paper. 

The axi-symmetric background should include both a Mestel dis c and an isotherma l 
halo. T he complete potential for such s tructu r e is well known (e.g. iBisnovatyi-Kogan 



( 19751 ). iMonet. Richstone fc Schechter I (198ll ). iToomre I ( 19821 )). and a self- similar 



derivation is given in (IHenriksen 



( 20121 )). In general the halo also possesses a non 
axi-symmetric component. However we ignore the non axi-symmetric halo here as it is 
likely to be a small effect at the disc. 



The axi-symmetric isothermal halo potential at the disc is of the form 



with constant, but it is removed in the locally rotating frame of the disc. The coefficient 

(r) 

gives the radial potential component of the initial self-similar spiral structure. We will 
find that it should be substantially larger than the non-axially symmetric potential 

The Mestel Disc potential takes the form (e.g. 

27rGS a 



Henriksen I ( I2012h ) 



8 



(InSr + ln(l + costf)) 



where we have set 5 = a/V and the axi-symmetric surface density is a a 
plane the total inertial potential becomes 

/ 2 7 rGE a 



(20) 

E a /(5r). In the 



$ = $ a (7r/2) + $ dr + $ dh = $ oa In 5r + ^ 



dr 



+ $£ + *J£)h <&• + **■■ (21) 



In the local disc reference frame however, we may focus on the non-axi-symmetric 
spiral potential (including the axi-symmetric part) since 



V 2 



2ttGS 



S 



a + $W 

1 do ' 



(22) 



This free-fall condition removes the influence of the undisturbed disc-halo potential on the 
spiral particles. The developing spiral structure perturbs the radial and angular structure 
of the disc through the action of $d r . 

The distribution function of the axi-symmetric particles is not our concern i n this paper 



but it is also approximated by an isothermal distribution in the rotating frame (IHenriksen 



( 20121 )). The total axi-symmetric surface density may also contain isothermal gas, for which 
a collisionless DF is not relevant. 



3. Potential-Disc Coupling 

The asymmetric potential found from the Laplace equation (f!8l) must obey the disc 
boundary condition (the axi-symmetric part is cancelled by — (l/r)do§ a \g =n /2) in the 
absence of gas as 

2irGa dr = —{d e q dr )\e = „ /2 , (23) 
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where 

a dr = I F(E' d ) dY R dY^ = 2k$ { £kJ-^. (24) 

This is equivalent to treating the disc as a volume density with a delta function 5(z) or 
5 (cos 6) /r in the Poisson equation. Integrating about z = 0or# = 7r/2 gives this boundary 
condition on the upper half space. 

This condition is generally difficult to satisfy for all k for a single value of m, because 
of the exponential dependence of a on the potential that follows from the last integral. 
Fortunately we can arrange to satisfy it everywhere by adding isothermal gas on the left of 
equation (1231) to obtain the augmented equation (1541 below. The gas is not described by 
the isothermal collisionless DF and so provides an additional degree of freedom. In principle 
its distribution is determined by the magnetohydrodynamic equations in the potential 
but we avoid this in this paper by simply using the augmented equation (|23l) as an equation 
for the gas surface density. 



Equation ( 1181) with the base ansatz of equation (11911 is resolved into two equations 
(e e/a for compactness) 

1 d / , ^ dT^Yf) x rrl / jn \ / * 2 2/2 



(sin0-f ) + T m (6)(k 2 m - m 2 (e 2 + —^-) + ime) 



sm9d9 K d6 m \ j \ m v sin 2 e < 
= 0, (25) 

d( 2 d( 
= 0, (26) 

where fc^ is the separation constant (positive or negative or zero) and ( = Vt/r. 
The base solution for the radial dependence takes the form 

K m (0 = expz(-^ + ^lnC)(A lm M AiM ( v / 3mC) 

+ A 2m W x ^(V3m()), (27) 



where M,W are Whittaker functions, A = y/3me/2, and /i = (a/(1 + ime) 2 + 4fc^)/2. 
One must remember that this asymmetric solut ion is to be added to an axially symmetric 
Mestel disc or disc-halo (e.g. iHenriksenl ( 20121 )). so that positive and negative values 



are acceptable. The amplitude must however be less than that of the axially symmetric 
potential. 

Our base solutions are not normal modes since they are not standing waves. However 
in the full Mestel disc without imposed boundaries they are the modes available to the 
disc. The various possible asymmetric potential base behaviours as a function of £ can be 
seen for typical values in figure([T|). We have applied the exponential phase factor e imK in 
equation ( JT9l) to the potential. One should recall (from equation ( JT2l) ) that we can write 
k = v + (, where v = <p + (e/a) lnr is the undisturbed logarithmic spiral for v constant and 
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£ = fit. We show only the real part and amplitude of the functions, but the imaginary part 
might also be relevant if the multiplicative factors are imaginary. 

We see from the figure that the 'M term' represents an oscillation growing in ( from 
zero amplitude. It is ultimately unstable and grows exponentially after ( « 8. The case that 
we illustrate has m = 2 and e/a = 3.5. A similar calculation with e/a = 2.5 reveals that 
this instability is similarly violent. It becomes more violent with significantly decreasing 
winding angle. 

However it is clear that for any t > this term is singular at the origin since ( will 
be infinite there. This singularity is evidently the formal source of the instability since 
the value at any fixed ( propagates to larger radius in time. However we must remember 
that we require C,Yr/V <C 1 for our potential solution to be valid. Thus there is only a 
limited range of ( accessible to our solutions and infinite ( is excluded unless Yr tends to 
zero at small radius. Generally there is an effective inner boundary at finite t, set by some 
maximum ( that depends on the maximum value of Yr/V. This also limits the growth in 
time to a maximum ( at a given radius unless Yr = everywhere. 

The potential term proportional to the W function' represents the development 
from an initial spiral of fixed v of finite amplitude. The growth is also oscillatory (in 
both amplitude and real value) but rather than growing exponentially, it is exponentially 
damped after passing through a strong localized amplification, which we identify with swing 
amplification. In contrast to the M term, we find that the amplification becomes rapidly 
weaker as the winding angle e/a declines. 

In this paper we focus on the W function behaviour since the origin of the initial arm 
may be left undetermined. The history that we study is then a transient swing amplification 
that ends in an exponential decay. The potential determines the arm surface density 
through equation (1231 . 

We choose the separation constant to be k m = as an initial condition that yields 
an initial spiral in the Whittaker W base mode. This is so because, provided that 
Re(2fi + 1) > 2 and u niquely when k m = 0, the term in the potential proportional to the W 
function goes to (e.g. lAbramowitzfc Stegun I ( 19721 )) 



A 2m T m (7r/2)e- imeln( ^ )/2 — -) + cstxzx e imeln W 2 ) e im », (28) 

r(/i-A + i/2) j 

as ( — > 0. Here the argument \/3m£ = z. The product of A 2m T m (7r/2) exp (— imeln (\/3)/2) 
with the ratio of the Gamma functions is a complex constant, wherein A is -y/3me/2 as 
above, and now /i = (1 + ime)/2 since k m = 0. This yields our initial condition of a 
pre-existing arm at ( = since the latter expression for the limiting W potential tends to 
this complex constant as ( — )■ 0. 

When k m = the solution for T m (8) is given in terms of the familiar associated 
Legendre functions as 

TJ9) = (C lrn PZ € {x) + C 2rn QT mt {x)), (29) 

where and denote the associated Legendre functions and C\ m , C 2m are complex 
constants. 
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Fig. 1. — On the upper left we show the real part of the asymmetric potential proportional 
to the whittaker M function, as a function of ( for v = 0, m = 2, and e/a = 3.5. The 
amplitude constant T 2 (7r/2)A 12 is omitted. The upper right curve is the modulus of this 
potential component. At lower left we show the real part of the potential proportional to 
the Whittaker W function for the same parameters (T 2 (7r/2)A 2 2 is omitted). At lower right 
is the amplitude of this potential component. The same graph as at upper left for a leading 
spiral with e/a = —3.5 diverges much more rapidly starting at ( ~ 0.7. In every case k m — 0. 
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The behaviour that interests us most is that for ( > 1, when the evolution of the 
initial spiral is marked . Usin g the asymptotic forms applicable when z = \/3m( > 1 



Abramowitzfe Stegun I (119721 ) the potential at the disc in this regime takes the form (we 



can absorb A 2m into the constants in T m j 

^ dr ~ T m (vr/2) exp im(v + C/2 + e In (/2)z x e - z/2 . (30) 

We see that the phase of the original spiral v is modified by the extra terms depending on 
C in the exponential plus any phase constant in T m (-7r/2). The amplitude is modulated in 
C by the function z x e~ z ^ 2 , which peaks at z = 2A or £ = e/a. El This last value should be 
> 0(1) for the peak to lie in the asymptotic range. The value at this peak is e A ( ln2A - 1 ) ) 
which can be extremely larg e. This is the gr owth that we identify with a non-linear version 
of the 'swing amplification' ( jToomrel ( 19821 )). 



The 'swinging' in the evolution can be seen by setting the change in the total phase 
of the potential to zero. This identifies the original spiral in the asymptotic region where 
it is evolving in time and space. Using the phase in equation ( 1301) we find from the zero 
differential 

^ = -f(l + -7). (31) 
2 a( 

We recall that dv = d<ft + (e/a)dr/r, and we have restored e to e/a. At fixed r, dv = dip and 
d( = Vdt/r, so that the previous equation implies 

dd> V e 1 . . . 

^ = -^l + - 7 - 32 
dt 2 aQ 

Hence an initially trailing arm (e/a > 0) is winding up in the trailing sense at each r with 
an approximate maximum rate at £ = e/a of —V. An initially leading arm (e/a < 0) swings 
from a winding forward rate to a winding backwards rate at the peak amplification where 
£ = | e/a |, and continues at large ( to reach the limiting trailing rate of — V/2 thereafter. 

The true shape of the spiral at any instant is extended in r at fixed t. From that 
perspective, d( = —((/r)dr and dv = d<p + (e/a)(dr /r), and hence equation (j3~Ij) yields 

T Tr = 2^-a ) - (33) 

This implies that an initially trailing wave (e/a > 0) is now leading at small r (large () but 
still trails at large r (small (). The transition is again at the peak value Qt = ( = e/a or 
at t = r(e/a)/V. An initially leading spiral (e/a < 0) is tightly wound at small r (large Q 
and is more open at large r (small (). We can not extrapolate strictly to very small ( since 
equation ( 1301) holds for large z = \/3m(, but either side of the peak value is permitted. 

This behaviour is indicated in figure (|2]), where the phase from equation ( 130]) is plotted. 
We see there that at small r we have indeed a leading spiral, which transits to a trailing 
spiral at ( = e/a, the swinging amplifier peak. An initially leading spiral unwinds slowly 
from the centre outwards as expected. 



We recall that e in the formulae denotes the initial spiral winding angle e/a. 
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Fig. 2. — The figure on the upper left shows a spiral with initial winding angle e/at = 3.5 and 
Vt = 1.75 in the xy plane. The arbitrary numerical value of the total phase is set equal to 
10 for the illustration. The figure on the right shows a leading spiral with the same winding 
angle at Vt = 5. The bottom left figure shows the same trailing spiral at the sequence of 
times Vt = 1.5, 1.75, 2.0. The spiral moves out and swings back as time advances. The 
bottom right figure shows the amplitude modulation function from equation fl30|) for m = 2 
and e = 3.5 as a function of £. 
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The independent self-similar evolution with a < and ( = —V\t\/r is complementary 
to the solution above. Considering the structure in radius, we see from equation (133]) that 
an initially trailing spiral {e/a > 0) trails everywhere, but most tightly at small radius. 
An initially leading spiral (e/a < 0) is trailing at small r and leading at large r, with the 
transition occurring again at ( = \e/a\. Thus the behaviours of the leading and trailing 
spirals are interchanged from the case with a > 0. This is also the case with the rates 
from equation f )32|) . The interchange of leading and trailing spirals with different sign of 
the winding angle e/a at small r, suggests the possibility of the spiral structur e s aris ing 
from waves passing through the centre of the system (e.g. iBinney fcTremainel ( 20081 ) .pp 
512-513). 

Such behaviour might correspond to the self-excited M term in the potential, if the 
spiral waves are reflected at spatial infinity where ( = and so the amplitude is zero 
there. It still requires Yr — y for the centre to lie within the range of our approximation. 
Moreover the swing amplification must operate on the leading and trailing waves. In fact 
the amplification seems to much more rapid for leading waves (see the caption of figure ([I])). 

We note that the inner spirals, either leading or trailing are distorted to the opposite 
winding at a radius r = Vt/(e/a). This radius is 10 kpc after 0.18 Gyr if the winding 
angle e/a = 3.5. The initial pitch angle is about 16°. This interval is just over one half 
of a rotation period at this radius, but the outer spiral is not destroyed by this time. The 
oscillating nature of the potential seen in figure ([3]) below corresponds to the winding of the 
arm in radius and the outwards propagation of the arm in time. When m = 2, any value of 
v plus v + 7r represent the two original arms. 

The bottom right curve in figure (T5]) shows the swing amplification amplitude when 
e/a = 3.5, and m = 2, as a function of (. The maximum is pronounced in this case, but it 
declines sharply with declining e/a and m. Thus when e/a = 1 (so that the pitch angle and 
the winding angle are equal at 45°), the maximum is ~ 1.52. At e/a = 2 (winding angle 
63.4° it is w 25.56. For m = 1 and e/a = 3.5 it is only « 11.37. 

In any case the amplitude is subject to multiplication by an arbitrary constant (the 
previous paragraph discusses relative variations), which should be sufficiently small that 
the asymmetric potential does not dominate the axi-symmetric potential. This can not be 
done indefinitely if the exponential M instability is really present. In such a case the growth 
must end by drastically rearranging the disc. This would involve some form of dissipation, 
probably due to gaseous shocks and star formation. 

An important property of the swinging amplification on an intial spiral as seen in the 
W function behaviour, is that it rises rapidly at £ = 1 and extends with non-negligible 
amplitude over a considerable range in (. Thus it amplifies over the interesting asymptotic 
range of the potential. 

We can obtain a plot of the spiral winding in ( for given v by taking appropriate cuts 
at constant values of Re(^dr)- This function is shown in figure ([3]) and is expressed in 
equation (1271) . The cuts for the W and the M components are shown separately in figure 
()3]). We ignore the constant coefficient T m (n/2). 



To use this figure one can choose a given z/, which then 'recurs' in ( (that is increasing 



2 4 6 8 10 12 2 4 6 8 10 12 14 



Fig. 3. — The figure shows on the left v as a function of ( when the real part of the potential 
is 1 for the W component. On the right we show the curves when the real part of the 
potential is set equal to 1 for the M component. We set m = 2 and e/a = 3.5 in all cases. 

time and decreasing radius). We see this recurrence several times corresponding to different 
windings. The value v + 7r then describes the winding of the second arm as a function 
of (. In the range v — > v + 2tt there are actually two possible pairs shown, but only one 
pair is relevant depending on the chosen v. The number of windings are finite for the W 
component because of the rapid dissipation, while any potential cut gives windings for the 
M component until the exponential explosion. Moving vertically in the figure is equivalent 
to moving in 0. The negative slope part of the arm reflects approximately the swinging 
of equation ( 13TI) . since near the maximum a cut at constant potential is nearly a cut at 
constant phase ( 1301) . 

It is interesting to note that the growth always begins first at small radius and then 
propagates outwards in every case. We find that by reducing the winding angle (increasing 
the pitch angle) the arms dissipate more rapidly in (. At e/a = 2.5 the arms begin to 
dissipate by ( =~ 4 rather than ^ 6 in the case illustrated in figure ([3]). The amplified arms 
are thus more concentrated in time and space. The amplification is reduced in amplitude 
for single armed trailing spirals with m = 1. 

The real part of the W component of the potential has much the same form for 
negative e/a, that is for leading spirals. However its amplitude is reduced by many orders 
of magnitude. The M component is exponentially unstable much earlier for leading spirals. 
This is illustrated in figure (jl]) for spirals with e/a = ±3.5 and m = 1. 
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Fig. 4. — The left figure shows the M potential component for a single armed spiral with 
e/a = 3.5. On the right is the same component of a leading single armed spiral with 
e/a — —3.5 The instability is much more rapid for the leading spiral. 



4. Disc Boundary Condition 



The boundary condition ( I2"3"j) can only be satisfied at discrete values of (, unless there 
is a collisional asymmetric gas component a gr present in the disc. Such values would 
approach self-consistent normal modes as closely as possible in the complete self-similar 
disc. Otherwise a modal analysis would have to be abandoned in favo ur of integrals ove r 
the base functions, preferably using spiral arm based coordinates (e.g. iHenriksenl ( 201l[ )). 



In this paper we avoid this analysis (although it is quite feasible) by adding gas that is 
after all present in real galactic discs. Assuming a gas contribution converts the condition 
(E3D to 

1 



2ttG(<7. 



dr)\e=n/2, 



(34) 



where the expression (1241) for Odr applies. The gas surface density must be compatible with 
the isothermality of the disc, including its spiral structure. We therefore assume that it 
takes the isothermal self-similar form 



a 



gr 



S flr(«,C)/(^)> 



(35) 



where it is convenient to take 5 = a/V. This converts equation ( 1231) into a normal mode 
analysis for the joint system of gas and disc. In general the gas dynamics is complex and 
involves the galactic magnetic field as well as pressure and gravity. However this complexity 
goes beyond the simple normal mode model. The present procedure requires the gas to 
follow an arm, to which potential the gas behaviour is in fact very sensitive. 
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We recall that the asymmetric gas density a gr in the condition (|34"|) can be negative, 
so long as its magnitude does not exceed that of the background axi-symmetric density a, 
The self-consistent condition fl34|) now becomes explicitly for a single mode m (x = cos 6) 



S%(*c, C) + l^Kl exp H *^*' )) = ^«(^,()| 1= o. (36) 

The potential and its derivative at the disc follow from equation f lT9"j) . where the real part 
of the potential is required. When ( is small we have the initial logarithmic sp iral and 71™ . 
may b e taken as 1 in equation (fT9~j) . This case has been discussed in reference iHenriksenl 
( 20121 ). but it is of less interest than the evolving arms to be considered here. In the 



regime of evolution £ > 1 the relevant solution is given by the real part of equation f[3"Uj) . as 
discussed above. 



Unlike the situation described in ( Henriksen I ( 20121 )). the solution for T m (x) given in 



( 129|) can take a simple form. In the earlier work we anticipated the possibility of using the 
term in Q'^ ne to cancel an infinity at the poles due to an axi-symmetric disc-halo potential. 
In the absence of such an infinity we may simply take C2m = and impose the desired 
symmetry according to P^ t (—x) = P^ t (x). 

The form of the spiral potential in our regime of interest for the evolving initially 
given structure is that of equation ( 1301 . Since k = v + (, the various functions are most 
instructively regarded as a function of v and (. We consider as our example the tightly 
wound two-armed case with m = 2 and e/a = 3.5. This gives the 'modal potential' 
explicitly at the disc as 

<M0) = lCul1 ^ lo zV 2 / 2 cos (2(1/ + C/2 + - In fl + 4>f (0) + 12 ), (37) 

where z = 2\/3£ and A = 7\/3/2. The phase of the associated Legendre function is 
0^ 2) (O) -0.8915 and its modulus is \P$ | w 4.449 x 10 5 . We have defined ip dr = *^V $ S- 
The value of the amplification factor z x e~ z ^ 2 at the maximum, is ~ 8.64 x 10 3 . 

We are free to choose the arbitrary phase 0i2 = — 4> P (0) to simplify the argument of 
the cosine in the expression for ipdr- This choice of phase renders the potential negative 
on the arm. If moreover as an example we take ipdr = 10~ 3 and the cosine function times 
the amplification factor to be 1, then the corresponding v(Q curves are shown at left in 
figure ([3]). This is because the amplification factor times the cosine function gives W in 
this regime, and figure (JSJ) is plotted for Re(W) = 1. Numerically the chosen value of 
ipdr, together with the value of modulus of the Legendre function, imply that we have set 
IC12IM? ^2.25 x 10- 9 . 

The straight line in figure ([3]) that passes through the point v = 10, ( ~ 6.458375, has 
the approximate equation 

C 7 ii 
2(u + I + - In C) » 21- - 4.7415 x 10" 4 , (38) 
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and we can use this to evaluate the potential and its derivative at this point on an evolving 
arm that was initially characterized by this value of v. Substituting relation (138]) into 
equation fl36l) then gives the transient surface density of the gas on the arm. 

To proceed in a less encumbered way we write condition fl36l) in the compressed form 

Sg(vX) = qd x Tpdrk x =o)-e~ fdr , (39) 



where 



(2) 



s - 



9 " ***$K$ 



a/V_ 
7( 



q = (40) 



Thus S g is the gas density essentially in units of the collisionless density. The parameter q 
is a ratio of the characteristic dynamic scale to a gravitational scale namely aLc/V, where 

the gravitational scale Lq = l/(^7r 2 GK^}). The ratio is similar in spirit to the Toomre Q 
value, if a/(2n) is identified with the epicyclic frequency k = \/2VL, V/(2%) replaces the 
velocity dispersion, and KV 2 replaces the surface density. 

The derivative of the potential takes the form at the disc 

flllMo = lCulldP ^ /&:| V e-/ 2 cos (2(1/ + C/2 + ] In C) + $ P + M, (41) 

where the modulus \dPffi/dx\o ~ 3.24 x 10 6 , and the phase ~ —0.9576 so 
that0^ + 0i2 w -0.062. 

We calculate the gas density on this arm [y = 10 initially) at ( = 6.4584 by inserting 
u(Q from equation (138!) into equation (|39|) . Some calculation using our various numbers 
gives S g ~ — 0.959g — 1 on this arm at the chosen (. Thus there is a deficiency in the gas 
density at this (. 

The variation of the gas density with ( that follows from equation (1391) is shown by 
the large curve on the left in figure (jSJ) for an arm with an initial v = 10. The calculation 
follows the line passing through the point (£ = 6.4584, v = 10) in figure ([3]), which is given 
by equation (|38|) ). This is the region where the initial arm is rapidly winding. 

We note that, starting from a large excess on the outside of the arm (small (-not fully 
shown in order to display the potential), the gas density drops to strongly negative values as 
one crosses the arm towards the inside(larger ( at fixed time). Ultimately, leaving the arm 
on the inside, the gas density rises to a modest positive maximum and finally finishes as a 
negative constant where the resonant function drops to zero. On the same dimensionless 
scale the smaller curve gives the variation of ipdr with (. The gas density minimum is 
slightly inward of the potential minimum in the arm. 
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It should be noted that the amplitude of the variation in the gas density is very 
sensitive to the amplitude of the potential. Increasing this by only one factor of ten renders 
the swing in the gas density quite unreasonable. Such strong dynamic behaviour would 
require a full treatment of shocks and magneto-gas-dynamics. Given m and the winding 
angle e/a, the only real free parameters are this amplitude and q, the analogue of the 
Toomre Q. Our examples should thus be representative of real behaviour. 

The right-hand side of the figure gives the angular dependence of the gas density at a 
fixed ( and radius (3.5 In ar/V = 10). The two peaks are separated by two slighter deeper 
troughs. Since the angle increases in the sense of rotation, one encounters the peaks on the 
outside of the arms. 



5. Conclusions 

This paper has studied finite spiral structure in isothermal (constant rotational 
velocity) discs. The collisionless particles that comprise these spiral 'arms' are described by 
an isothermal distribution function in the locally rotating frame. The particles, the arm, and 
indeed the potential, therefore comove. 

The presence of gas is a necessary feature. We calculate its distribution based solely 
on the requirement that the thin disc self-consistent condition (|39|) be satisfied. The gas 
peaks and valleys are associated respectively with the potential peaks and valleys, so in this 
approximation it is also comoving. The gas density is dramatically higher on the outside of 
an arm. 

These structures are not steady and we found the time dependent potential that 
describes their evolution. For each 'mode' (we use this expression for the potential that 
satisfies equation (1391) for each m), this potential possesses two growth behaviours as 
shown in figure (Q. In the Whittaker W swing amplification studied mainly in this paper, 
the evolution of an initial arm is followed in time and space. This arm is first amplified 
as it 'swings' in pitch angle and finally dissipates exponentially. The 'swinging' and the 
amplification function are illustrated in figures (T5]) and ([3]). We use a winding angle of 74° 
(that is a pitch angle of 16°) for a two-armed spiral. This is tight winding. The argument 
is readily extended to other winding angles and normal modes (m). 

The disturbance in the gas density that follows by requiring the self-consistent 
condition (1391) to be satisfied everywhere on the disc, is shown in figure© for an initial arm 
undergoing swing amplification. The gas density is increased on the outside of the arm, 
declines strongly on the potential arm, and rises slightly inside of the arm. It finally attains 
a constant negative value inside the arm as the potential declines exponentially. At a fixed 
radius the same sequence is encountered in time as the arm crosses that radius. Such a 
sequence must be restarted in some undefined fashion. The variation in angle at a fixed 
value of the self-similar variable zeta is also shown in that figure. 

The amplitude of the variation in the gas density becomes exaggerated relative to 
any reasonable background, if the ratio of the asymmetric potential amplitude is greater 
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Fig. 5. — On the left the gas density that satisfies the boundary condition ( )39|) for v = 10, 
m = 2 and e/a — 3.5 is shown as a function of £ by the large curve. The small curve shows 
ipdr- In each case we start with the rapid winding of the arm at £ = 6.458. As £ increases 
with decreasing r at a fixed time, the plot shows the gas density falling steeply to a minimum 
from an excess (not shown) on the outside of the arm, then rising to modest positive values 
as we leave the arm on the inside, and finally attaining a constant negative value inside the 
arm. The figure at the right shows the variation of the gas density in increasing angle (that 
is with increasing v at fixed radius) at a fixed ( = 6.458. In each case q = 1.5. The peaks 
are on the outside of the arms. 
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than about 10~ 3 times the radial potential amplitude (essentially $^). In reality, stronger 
disturbances would require solving for the actual magneto-gas-dynamics in the asymmetric 
potential. 

The 'self-excited' instability of figures ([H |H) seems t o duplicate behaviour that which 



occurs in simulations of an isolated disc f Sellwood~l ( 201lh ). Its origin in a central singularity 



renders its physicality suspect. However it may indicate the presence of a genuine instability, 
which is due to the interchange of leading and trailing waves across the centre of the system. 
All that one can say definitely is that the amplification of both leading and trailing waves 
originates at the centre. 

Taking it at face value, the exponentially growing disturbance must be stabilized in 
some dissipative fashion, but we can not follow this development in general beyond a 
maximum value of £ = Vt/r. For completeness we see from figure (J4j) that the instability 
is slightly less rapid for a one-armed spiral than that of a two-armed spiral of the same 
winding angle. However it is much faster for the one-armed leading spiral of the same 
winding angle. 

In the case of a larger number of arms, such as m = 3, one finds (not shown here) the 
self-excitation to be similarly strong for the same winding angle. The swing amplification 
of an initial arm is of much larger amplitude however, unless the winding angle is reduced. 

On the whole this analysis suggests that galactic spirality is transient and either 
continually self-excited after each cycle or re-excited by external influences. It has a 
dramatic effect on the gas distribution in the disc, which is in turn essential to the 
continuous spirality. 
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